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Abstract 

We report the prediction of temperature-dependent diffusion coefficients of interstitial liydrogen, 
deuterium, and tritium atoms in a-Ti using transition state tlieory. The microscopic parameters in 
the pre-factor and activation energy of the impurity diffusion coefficients are obtained from first- 
principles total energy and phonon calculations including the full coupling between the vibrational 
modes of the diffusing atom with the host lattice. The dual occupancy case of impurity atom in 
the hep matrix are considered, and four diffusion paths are combined to obtain the final diffusion 
coefficients. The calculated diffusion parameters show good agreement with experiments. Our 
numerical results indicate that the diffusions of deuterium and tritium atoms are slower than that 
of the hydrogen atom at temperatures above 425 K and 390 K, respectively. 
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Most metals and alloys are susceptible to hydrogen damage, including phase transforma- 
tions, hydrogen embrittlement, and corrosion. Diffusion of hydrogen in metals provides an 
intriguing and challenging intellectual problem in addition to the practical impact. Taking 
into consideration the irnportance of diffusion process, many experimental measurements 
have been performed [l y]. Also, there have been many theoretical efforts to gain diffu- 
sion coefficients for multiple materials and the corresponding diffusion mechanisms using 
fundamental electronic or atomistic approaches |7-10|. In fact, although the diffusion ex- 
periments can determine the overall diffusion coefficients, they cannot generally determine 
the microscopic physical processes involved in the diffusion steps which are quite important 
for basic understanding and practical applications. Therefore, first-principles calculations, 
which have been widely used in the study of solid-state diffusion, can help to track the 
microscopic diffusion processes and provide specific quantitative values, such as formation 
and migration energies, involved in the diffusion process. 

Metal titanium (Ti), which combines high strength and low density, is a suitable material 
for the application in aircraft construction and aerospace engineering, as well as in the 
chemical industry. In most cases, titanium and its alloys display excellent resistance to 
damage. However, due to the strong affinity between titanium and hydrogen that results in 
a severe deterioration in the mechanical properties, considerable restrictions exist in their 
applicability under hydrogen-containing environments, where the hydrogen can be supplied 
by a number of sources, including water vapor, picking acids, and hydrocarbons. Despite 
the extensive work on diffusion, there is a serious lack of reliable experimental diffusion data 
on titanium and many other important systems. It is well known that metals of the IV 
group exhibit phase transition from a-phase (hep structure) to /3-phase (bcc structure) at 
the reaction temperature, which is 1155 K for Ti. This relatively low transition temperature 
for Ti makes diffusion experiments difficult to be solely carried within the a-phase. To date, 
the experimental diffusion coefficients of hydrogen in a-Ti are scattered over many orders of 
magnitude [iml], while to our knowledge, systematic ab initio studies of hydrogen diffusion 
in Ti is still lacking in the literature. 

The purpose of the present work is to comprehensively investigate the atomic diffusion 
mechanism of hydrogen and its isotopes (deuterium and tritium) in a-Ti, and to obtain 
the specific values of the corresponding energy barriers, vibrational free energies, and diffu- 
sion coefficients from first principles. In the most general form, the diffusion coefficient is 
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expressed as 



D = Doe 



-Q/kT 



(1) 



where Q is the activation energy, Dq is a pre-factor, and k is the Boltzmann constant. 
According to the transition state theory (TST) ll|, |l2|, the jump rate is written as 
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where Zts and Zq are partition functions for the transition state and the ground state, 
respectively, h is the Planck's constant, and AHm is the enthalpy difference of migration 
between the transition state and the ground state including the thermal electronic contri- 
bution. Within the framework of the harmonic approximation, the quantum mechanical 
partition functions in Eq. (2) can be expressed as 
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where uf^ and are the vibrational frequencies at the transition state and the ground 
state, respectively. Using the expression of the phonon free energy 
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the jump rate can be simply expressed as 

kT 
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where the zero-point energy is included in the Fyn, term. 

As depicted in Fig. 1, there are two stable interstitial positions, i.e., the tetrahedral 
(T) site and octahedral (O) site, respectively. In the present work we focus on the dual 
occupancy case. Thus, four jump frequencies are to be considered. In the hep-structured 
a-Ti, there are two kinds of impurity jumps, viz. one perpendicular to c axis and the other 



along the c axis, resulting in two diffusion coefficients which can be expressed as 
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FIG. 1: The arrangement of interstices in the hep lattice (left panel) and the energy profiles for a 
hydrogen atom to diffuse along (a) 0^0, (b) T^T, and (c) O— )-T and T— t'O paths. 

where u^o, wot, i^oo; and wtt are respectively the jump rates along the four paths shown 
in Fig. 1. 



The density functional theory 



DFT) 



calculations are carried out using the Vienna ab- 
initio simulation package (VASP) 3, 15| with the projector-augmented- wave (PAW) poten- 
tial method [16,]. The cutoff energy for the plane- wave basis set is 450 eV. The exchange and 



correlation effects are described by t he g eneralized gradient approximation (GGA) in the 
Perdew-Burke-Ernzerhof (PBE) form 17| • We employ a 3 x 3 x 2 a-Ti supercell containing 
36 host atoms to simulate hydrogen migration in the a-Ti matrix. To check the convergence 
of the formation and migration energies, we have also considered a 4 x 4 x 3 a-Ti supercell 
containing 96 host atoms and one hydrogen atom. The integration over the Brillouin zone 
is carried out on 3 x 3 x 3 k-point meshes generated using the Monkhorst-Pack ISj method, 
which proves to be sufficient for energy convergence of less than 1.0x10"^ eV per atom. 
During the supercell calculations, the shape and size of the supercell are fixed while all the 
ions are free to relax until the forces on them are less than 0.01 eV A"^. The phonon spectra 
are carried out using the density functional perturbation theory. 

The theoretical equilibrium lattice parameters a and c/a for a-Ti are 2.939 A and 1.583 A, 
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FIG. 2: Calculated phonon dispersions of a 36-atom a-Ti supercell with hydrogen, deuterium, and 
tritium atoms in the stable octahedral site (Oct), the transition states along the 0^0 (TS-OO), 
O— >T (TS-OT), and T^T (TT-TS) paths, and the metastable tetrahedral site (Tet), respectively. 
The modes with imaginary frequencies at the transition states correspond to the motion of the 
hydrogen, deuterim, and tritium atoms across the barriers. 

respectively, in accordance with the experimental data of 2.950-2.957 A and 1.585-1.587 
2l|. By comparing the formation energies of a hydrogen defect in 0-site and T-site, we find 
that the 0-site is more favorable in energy, lower than the T-site with a 36-atom supercell by 
0.080 eV. The 96-atom supercell gives the essentially same result with an energy difference 
of 0.082 eV. For the calculation of migration energies and phonon frequencies, an important 
issue concerns the location of the saddle-point position for hydrogen atom. Here, each 
saddle-point structure and the associated minimum-energy pathway (MEP) were calculated 
by employing the climbing image nudged elastic band (CINEB) method [22|]. As displayed 
in Fig. l(a)-(c), when hydrogen atom diffuses from the stable site to the transition site (TS), 
the energy profile can be well described by a sinusoidal curve, as originally suggested by Wert 



and Zener 



23|. In the case of O— )-0 path and T— )-T path, the energy is found to display a 
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TABLE I: Calculated migration enthalpy difference of atomic hydrogen diffusion via O— >0 

path, T->T path, 0-)-T path, and T^O path for 3 x 3 x 2 and 4x4x3 hep Ti superceU models, 
respectively. 
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single maximum, corresponding to a saddle point at the high-symmetry position located half 
way between neighboring sites. When hydrogen atom diffuses along the O— >T path or T— >0 
path, the transition state locates close to the T site, and the distance between the saddle 
point and T-site (0-site) is 0.852 A (1.001 A). In Table I, we listed the migration energies 
of the four diffusion paths using 36- and 96-atom models, respectively. By comparison of 
the two models, one can see that a 36-atom superceU is already convergent. Actually, the 
difference in migration energy of these two models are 0.04 eV at most for O— >0 path, 
indicating that a 36-atom superceU can be reliably used in the simulations of the single 
hydrogen difusion energetics. Thus, in the follows we adopt a 36-atom superceU to perform 
the phonon dispersion calculations. 

The phonon dispersions of the 36-atom superceU with a hydrogen atom at the 0-site, T- 
site, and three TS sites are shown in Fig. 2. When a hydrogen atom occupies an octahedral 
site, its phonons show a threefold-degenerate dispersionless branch at 30.5 THz. At the 
T-site, the frequencies of the hydrogen atom show a similar character with respect to those 
in the 0-site, and the threefold-degenerate branch moves to 42.4 THz. In both cases, the 
frequencies are positive, indicating a true minimum in the energy. When the hydrogen 
atom moves to a transition state, its phonon branches, as shown in Fig. 2, arc split into a 
degenerate doublet, which shift upwards to some extent, and a non-degenerate singlet, which 
shifts downwards so as to be imaginary, implying the spontaneous diffusion instability for the 
defective hydrogen in the transition state. By definition, the transition state is characterized 
by the occurrence of one negative eigenvalue in the dynamical matrix. Here the doublet 
modes at positive frequency correspond to the vibrations of the hydrogen atom perpendicular 
to the diffusion path, while the singlet mode at imaginary frequency corresponds to the 
motion of the hydrogen atom along its diffusion path. 
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TABLE II: Calculated diffusion pre-factors {Dqj_ and Do\\) activation energy {Q± and Qy) 
for atomic hydrogen, deuterium, and tritium in a-Ti. The temperatures represent the ranges over 
which diffusion coefficients were fit to extract corresponding pre-factors and activation energies. 
For comparison, experimental results are also listed. 



Method 


Do± (mVs) 


^011 (mVs) 


(kJ/mol) 


Q|| (kJ/mol) 


T(K) 


Hydrogen 


2.506x10"^ 


1.275x10-6 


51.90 


49.89 


280-1000 


Deuterium 


1.587x10-6 


8.630x10-'^ 


50.30 


48.50 


280-1000 


Tritium 


1.219x10-6 


6.781x10-'^ 


49.59 


47.85 


280-1000 


Expt." 


1.672x10-6 




51.80 




773-1097 


Expt.^ 


3.057x10-6 


r—i 


62.46 




884-1102 






" Reference [l]. 


^ Reference [3] 







When the hydrogen defect locates at the transition state along the T— J-T path, we note 
that the hydrogen-related branch has the highest frequency of 59.3 THz compared to 49.2 
THz in O— i-O path and 48.9 THz in O— J-T path. In fact, the transition site along the T— J-T 
path locates at the geometric center of triangle formed by its three nearest Ti atoms, thus 
the interstitial space is much more confined with respect to the other two transition states. 
Besides, there is a phonon branch between 9.0 and 10.0 THz in the transition state along 
the T— T-T path, which is obviously off from the bulk of Ti-related phonon branches. This 
branch is related to motions of Ti atoms, which couple with motions of the H impurity. 
As we expect, the impurity atom related frequencies decrease with increasing mass from 
hydrogen to its isotopes (cf. Fig. 2). 

Integration of the phonon dispersions over the entire Brillouin zone allows to compute the 
temperature-dependent enthalpy, entropy, and free energy, which can be used to evaluate 
the diffusion coefficient according to Eq. (1). In Fig. 3 we showed the Arrhenius plot 
of the computed diffusion coefficients together with experimental values in Refs. Illlsi for 
comparison. The diffusion coefficient curves for D± and D\\ are almost indistinguishable 
below 500 K. With increasing temperature, the value of D± becomes somewhat higher than 
D\\. Both of these two diffusion coefficients are in agreement with the experimental data by 

n 

Wasilewski et al. [1] and Papazoglou et al. [2;] (cf. Fig. 3). By linear fitting of the diffusion 
coefficient curves, we obtained the diffusion pre-factors (-Do± and Z^oll) and activation energies 



7 



800 



Temperature (°C) 
300 200 100 



10" 
10"' 
10"' 



(/3 



10 



10 



10 



■11 



10 



10 



■14 



15 



-16 




• Wasilewski et at. 

■ Papazoglou et al. 

■ Johnson et al. 



10 

1.0 1.5 

1000/T (K"') 

FIG. 3: (Color online) Calculated and experimental diffusion coefficients of atomic hydrogen, 
deuterium, and tritium in a-Ti. The Ttr=1155 K stands for the phase transition temperature of 
metal Ti from a-phase to /3-phase. 

{Q\_ and Q\\) according to the Arrhenius plot of Eq. (1), i.e., ln(D) = In(Do) — Q/kT. In 
Table 2 we listed our fitting results of pre-factors and activation energies. A linear fit between 
280 K and 1000 K of the calculated data for hydrogen gives (5^=51.9 kJ/mol and Q||=49.89 
cJ/mol, in agreement with experimental data ranging from 51.80 kJ/mol to 62.46 kJ/mol 
l|,|2|. The pre-factors are Do±=2.506 xlO^^ m^/s and Do||=l-275xlO~^ m^/s, respectively. 
One can see that Dq\\ is nearly the half of -Do±- The value of -Do± is well located in the 
experimental range from 1.672x10"^ m^/s to 3.057x10"^ m^/s. 
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FIG. 4: (Color online) The ratios of the diffusivity of deuterium and tritium to hydrogen atoms. 
(_L) and (||) represent the impurity diffusions perpendicular to and parallel to c axis, respectively. 

The temperature-dependent diffusion coefficients for the isotopes deuterium and tritium 
in a-titanium can be calculated based on the electronic structure results of hydrogen in 
titanium. One would expect a lower diffusivity for the heavier isotopes, since the impurity- 
related frequencies decrease with increasing mass (Fig. 2). In Fig. 4 we depicted the 
ratio of the diffusivity of deuterium and tritium to hydrogen. One can see clearly that 
this is only the case at temperatures above 425 K for deuterium and 390 K for tritium. 
At lower temperatures, the ordering is reversed. The reason is that at low temperatures, 
the zero-point energy effect is more pronounced [3], which increases the jump rate through 
an exponential dependence, since the energy difference between the transition and ground 
states shows a reduction trend from hydrogen to its isotopes. As the temperature increases. 
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the zero-point energy effect is offset due to the lower frequencies of isotopes. At 425 K, 
these two effects achieve a balance for deuterium. Also, a crossover temperature of 390 K is 
predicted for tritium, as depicted in Fig. 4. 

In summary, we have systematically studied the temperature-dependent diffusion coeffi- 
cients of atomic hydrogen and its isotopes in a-Ti using transition state theory with accurate 
first-principle total energy and phonon calculations. The quantitative agreement between 
the calculated and experimental diffusion coefficients has been achieved, which demonstrates 
the validity and practicability of the present theoretical method. Our calculated ratios of 
the diffusivity of isotopes and hydrogen provide a good reference for future experimental 
measurements on a-Ti. 

This work was supported by NSFC under Grant No. 51071032 and by CAEP Foundations 
for Development of Science and Technology under Grant No. 2011A0301016. 
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